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The objective of this task is to calculatfe damage evolution 
at the critical location of selected components of reusable 
space propulsion systems as a function of their operating 
environment (s) and to relate the evolution of damage at the 
critical locations to the performance variables of the 
propulsion system using advanced constitutive and damage 
models . The models so developed will be oriented toward use 
in an advanced diagnostic/prognostic health monitoring 
svstem for reusable propulsion systems. Further, the models 
developed must be amenable to a probabilistic interpretation 
in the future. The intent is to build on the state of the 
art of previous work in these fields whenever possible 
rather than to develop new models. The components to be 
modelled include, a Lox Post (injector element) , a thrust 
chamber, and a turbine blade. The study conducted on the Lox 
Post is reported here. 


The loads that contribute to the main injector lox post 
damage evolution are temperature, static and dynamic 
pressure loads, and vibration loads. The temperature load on 
the post is controlled by engine system variables such as 
preburner mixture ratio, hot cas temperature and flow rate, 
Lox temperature and flew rate, local hot gas and coolant 
film coefficients, and any other local geometric and 
material property parameters . The pressure loads are a 
function of fluid drag and turbulence and static pressure 
differential. In this particular application damage at a row 
12 Lox post thread root location is computed and related to 
system variables. It is apparent that generation of a link 
between performance variables and damage will span several 
disciplines such as engine performance, fluid mechanics, 
heat -transfer , structural analysis, and material science. 


In general, for any component only a subset of the system 
variables that influence damage are directly measured in an 
actual flight. Thus it becomes necessary to rely on a 
numerical engine model that will compute the system 
variables that influence damage directly from a set of 
measured more primitive (independent) and/or local system 
variables that control engine health, thrust and 
performance. The most frequently used form of the model is 
the influence coefficient model which derives its origin 
from the more computationally intensive engine balance model 
and test measurements. Typically, the influence coefficient 
model is generated such that it accurately portrays engine 
performance in the plus or minus 3 sigma range of engine 
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performance . One such influence coefficient model 
part of this task. 


is used as 


The development of the system variable to structural 
response link require a large number of one time heat 
transfer and structural analysis runs covering the entire 
start transient, steady state, and cut off transient. 

Several par am etric heat transfer analyses were run to set up 
influence coefficient models for different regions of the 
Lox Post that help to quickly compute the entire temperature 
field in the Lox Post over the full duty cycle given the 
local system variables. Using heat transfer influence 
coefficient models several parametric nonlinear structural 
analyses over two duty cycles were run to establish the link 
between system variables to temperature, stress levels, 
strain levels and, multiaxial factors at the critical 
locations of the hysteresis loop. A very refined 
axisymmetric model was used for this purpose. 

Since the temperature was low enough for this component, the 
standard rate independant plasticity model was considered 
adequate for this application. 

The dynamic response due to random pressure and mechanical 
vibration was computed using a detailed beam finite element 
model. The pressure P.S.D. shape and power, correlation 
length and their relation to local mass flow rate and fluid 
density were obtained from fluid dynamics unit. These system 
variables are linked to other system variables through 
engine system models. The r .m. s . stresses due to mechanical 
vibration (harmonics and random) were obtained using a 
multibase excitation analysis. Thus the r.m.s. high 
frequency stress response due to dynamic excitation is 
linked to system variables. 


Thus at the end of above steps, the link between system 
variables to response variables that control damage 
evolution has been established. The damage computation 
module has been written with enough modularity to compute 
damage due to linear as well as nonlinear accumulation 
algorithms using the damage curve approach. Given the system 
variables, power level, and their duration the module 
computes the damage accumulated and if necessary computes 
the remaining flight life assuming that the last flight 
profile will be repeated in future flights. Several further 
options are available for the user to experiment with 
different stress amplitude distributions, order effect, 
hydrogen embrittlement and mutiaxiality factors, etc. 


Future efforts involve the examination of system variables 
on life prediction of Main Combustion Chamber Liner and 
Turbine Blades. The main combustion chamber liner effort 
involves the implementation of thermal rachetting and crack 
growth prediction models and turbine blade effort involves 
implementation of life prediction methodologies for 
anisotropic superalloys. 


LIFE PREDICTION TASK OBJECTIVES 


• DEVELOP LINK FROM PERFORMANCE VARIABLES TO DAMAGE 

• USE ADVANCED CONSTITUTIVE AND DAMAGE MODELS 

• MODELS USED MUST BE AMENABLE FOR EXTENSION 10 PROBABILISTIC DOMAIN 

• DEMONSTRATE METHODOLOGY FOR I.OX POST, MAIN COMBUSTION CHAMBER LINER, TURBINE BLADE 


LIFE PREDICTION LINK TO HEALTH MONITORING SYSTEMS 


I I 
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I I 
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I I 

I NUMERICAL ENGINE I • 

I MODEL I 
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I RESPONSE HISTORY I 
I COMPUTATION MODEL I 
I AT CRITICAL LOCATIONS I 
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! DAMAGE CALCULATION I 
I MODEL I 
I I 
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DAMAGi: CALCULATION USING ADVANCED MODELS 


221 




LOADS ON MAIN lN.ILClOR HI.CMENT 

ROW 12 
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Pressure Loads 


Critical Thread Area 


• DOMINANT LOADS 

e TCMPERATURE DIFFERENTIAL 

• FLUID LOADS (DRAG AND DYNAMIC PRESSURE r I.IICUIATION) 

• LESS DOMINANT LOADS 

• MECHANICAL V1DRA1ION 

• CRITICAL AREA LOOKED AT IS THE TlIRfAl) IOCA1IOH 


AN APPROXIMATE ENGINE PERFORMANCE MODEL 
INFLUENCE COEIIICIENT MODEL 


• EXTRACTED FROM NONLINEAR BALANCE MODI L 

• POLYNOMIAL REGRESSION FIT 

• STRONGLY CORRELATED WITH POWER LEVEL 

INFLUENCE COEFFICIENT - IC,j(PL) = C'll -I- C\PL + C>PL~ + C,\PL 
dependent variable mean value - V,(PI‘) = "ft + a ' PL + a - PL + a ' PL 

. magnitude evaluation or DEPENDENT variaiiies 
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>01! UM’dATtlti 


AN EXAMPLE OF INDEPENDENT AND DEPENDENT VARIABLE SET 


INDEPENDENT VALUES 


1. FACILITY MIXTURE RATIO 6.0260 

2. FUEL INL TOTAL PR (PSIA) JO. 0000 

3. OX ID INL TOTAL PR (PSIA) 100.0000 

NI FUEL INL TEMP (DEG R) 37.0000 

5. OXID INL TEMP (DEG R) 16'l.0000 

6. HPFTP TURB NOZZ AREA 10.8000 

7. HPOTP TURB EFF MULT 1.03)0 

8. HPOTP TURB NOZZ AREA 2.8960 

9. THRUST CHAMB C* MULT 1 0053 

10. HPFTP TURB EFF MUI.1 .‘JZ'iO 

11. POWER LEVEL I.O'lOO 

12. ETC. ... 


DEPENDENT VALUES 


1. HPOTP SPEED (RPH) 

28335.8600 

2. IIPF1P SPEED (RPM) 

35660.2700 

3. HPOTP PUMP DIS PR (PSIA) 

4321.3460 

<1. IIPF1P PUMP DIS PR (PSIA) 

6521.7120 

5. OPB CHAMBER PR (PSIA) 

5308.2630 

6. FPB CHAMBER PR (PSIA) 

5213. <670 

7. ENGINE OXID FLOWRATE 

930.5878 

8. ENGINE FUEL FL0WRA1E 

15*1.2625 

9. ENGINE THRUST 

490784.4000 

10. OXID PRESSHRAN1 F/R 

1.6580 

11. HILL PRESSURAN1 F/R 

.7280 

12. OPOV POSITION 

.6706 

13. FPOV POSITION 

.8001 

1H. MCC INJECTOR END PR 

3126.2170 

IS. IIPOP INLET PR (PSIA) 

386.2080 

16. IIPFP INLET PR (PSIA) 

206.8424 

17. PBP D 1 SCH PR (PSIA) 

7389.5130 

18. IIPOP INLET TEMP 

169.7184 

19. IIPOP DISCH TEMP 

192.0415 

20. IIPFP DISCH 1EHP 

98.1692 

21. PBP DISCH TEMP 

205.3588 

22. IIPFP INLET TEMP 

42.3058 

23. LP01P SPEED 

5166.0650 

2'l. I.PFTP SPEED 

15130,7500 

25. IIP01 DISCH TEMP A 

1285.8410 

26. IIP01 DISCH TEMP D 

1285.8410 

27. IIPM DISCH 1EMP A 

1735.9330 

28. ETC. ... 



THERMAL TRANSIENT ANALYSIS USING A REFINED THERMAL MODEL 


• ANSYS MODELING/SINDA THERMAL ANALYSIS • FOUR DISTINCT ZONES 
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ORIGINAL PAGE IS 
OF POOR Ql!ALl T Y 


AN APPROXIMATE THERMAL ANALY7ER MODEL 
QUICK COMPU1 AT IONS ON TEMPERATURE FIELD 


• INDEPENDENT VARIABLES 

• HOT GAS TEMPERATURE AND FLOW RATE • HOT GAS AND LOX h FACTOR 

• LOX TEMPERATURE AND FLOW RATE • HAYNES 188 K FACTOR 

• LOX PRESSURE 

• INFLUENCE COEFFICIENT MODEL 

• COEFFICIENTS OBTAINED FROM LARGE NUMBER Of PARAMETRIC RUNS 

Max.Twg — Ref .Mar-Ting IV — Ref.V \V — Ref.V} 2 
Ref.h4ax.Twg 1 Ref.V + ^ Ref.V 

Min.Twg - Ref.Min.Twg _ f \V - Ref.V ] f \V - Ref.V] 2 

Ref.h4in.Twg 1 | Ref.V j 2 Ref.V 


• DIFFERENT INFLUENCE COEFFICIENTS FOR FOUR /ONES 


STATIC NONLINEAR ANALYSIS FOR TWO DUTY CYCLES USING A REFINED MODEL 


• ANSYS COMPUTER CODE USED 

• SEVERAL PARAMETRIC ANALYSES CORRESPONDING TO PERTURBATIONS OF SYSTEM VARIABLES 


• TWO DUTY CYCLE TEMPERATURE CYCLING 



A REFINED MODEL A TYPICAL AXIAL o - € RESPONSE 
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ORfGIfMAL PAGE iS 

OF POOR QUALITY 



Psd (accel /Hz) 


APPROXIMATE STATIC ANALYSIS MODEL UNKING RESPONSE TO SYSTEM VARIABLES 

l.CF DAMAGE DRIVERS 


• RESPONSE VARIABLES 

• EFFECTIVE STRESS 

• TOTAL EFFECTIVE STRAIN 

• PLASTIC STRAIN 

• TEMPERATURE 

• MULT I AXIAL I TY FACIOR 

• TRACKED AT THREE POINTS OE HYSTERESIS I OOP 10 CALCULATE LCF DAMAGE 

• A TYPICAL INFLUENCE MODEL 

e to1 = bo -f Y, at. t he j ih point, of the loop 

3 i=l ,9 


DYNAMIC ANALYSIS OF LOX POST 
IICF DAMAGE DRIVERS 
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UNKING RMS RESPONSE 10 SYSTEM VARIABLES 


SCALING RELATION 


"Iwsin.n) = Pj f (^) 


. \3 


TOTAL 


0RMS ~ 'Ji a 'lt\lsY + [ a tiffs) 


• SINCE DYNAMIC RESPONSE IS DOMINATED BY I I RSI BENDING MODE OF HIE POST (1800IIZ), 
EXPECTED FREQUENCY HAS NO SENSITIVITY 10 SYSTEM VARIABLES 


DAMAGE ACCUMULATION 


• LOW AND HIGH CYCLE DAMAGE ACCUMULATION; LINEAR OR NONLINEAR (DAMAGE CURVE) 

• LCF STRESS-STRAIN STATE DEPENDENT ON SPECIFIC MISSION INPUTS 

• HCF DYNAMIC STRESS DEPENDENT ON 1IMF. AT DIFFERENT POWER LEVELS DURING MISSION 

• TEMPERATURE DEPENDENT FATIGUE CURVE 

• ADJUSTMENT OPTIONS TO INCREASE LCF STRAIN RANGE DUE TO HCF "MAX" AMPLITUDE 

• ADJUSTMENTS TO FATIGUE CURVE FOR MUM I AXIAL I TY AND HYDROGEN EFFECTS 


HCF ACCUMULATION 


• DYNAMIC STRESS AMPLITUDE DISTRIBUTION COULD BF. RAYLEIGH OR USER DEFINED 

• RANDOMIZED AMPLITUDE SELECTION OR BIN 1N1EC.RA110N ME 1 HOD TOR CALCULATING DAMAGE 

• SHAKE DOWN OF MEAN DUE TO PEAK STRESS EXCEEDING YIEI D 

• NONLINEAR MEAN STRESS CORRECTION FOR POSIIIVC MEAN STRESS 
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Mr or P ' (ksl) 


ADDITIONAL ELEMENTS Of THE DAMAGE MODEL 


• DETERMINATION OF CONSISTENT CONSTITUTIVE AND FATIGUE CURVE 

ASSUMPTIONS: RAMDERG-OSGOOD FORM OF CYCLIC STRESS-STRAIN CURVE. HANSON-COFFIN TYPE 

OF FATIGUE CURVE 


E 

1,0 F ' 


[a,, - (M‘ + c K )j\r [|^ + c (<ff - W * on;) 


min 


• EFFECT OF MAXIMUM RANDOM IICF LOAD ON I.CI STRAIN RANGE THROUGH NOTCH PLASTICITY 
ASSUMPTIONS: NEUBER'S RULE APPLIES AT POINT 2 OF THE HYSTERESIS LOOP 


O a,lj {t a ' b 


tf") = 


{a, + 3(j ru ‘" ,om )' 

E 


• EFFECT OF THE MULT IAX1AI.1TY FACIOR (ME > ON LIFE 
MODIFIED FATIGUE CURVE CONSTANTS: 



• EFFECT OF II, ENVIRONMENT THROUGH DUCTILITY LOSS 

MODIFIED FATIGUE CURVE CONSTANT ACCORDING TO UNIVERSAL SLOPES EOUATION: 



ADDITIONAL ELEMENTS OF THE DAMAGE MODEL (CONTINUED) 


• CUMULATIVE MULTI AXIAL AND ll 2 EFFECT IS DOMINATED BY MULTIAXIALITY NEAR THE END OF 
THE DUTY CYCLE 


Variation of Parameters durinf Startup Transient 
(second duty cycle, at the root of the critical thread) 

son 
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Variation of Parameters durinf Shutdown Transient 
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Temperature (F*) 




